Measurement of the Cavity-Loaded Quality Factor in Superconducting Radio-Frequency Systems 
with Mismatched Source Impedance 
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The accurate measurement of parameters such as the cavity-loaded quality factor (QL) and half bandwidth 
(fo.5) is essential for monitoring the performance of superconducting radio-frequency (SRF) cavities. However, 
the conventional “field decay method” employed to calibrate these values requires the cavity to satisfy a “zero- 
input” condition. This can be challenging when the source impedance is mismatched and produce nonzero 
forward signals (V;) that significantly affect the measurement accuracy. To address this limitation, we developed 
a modified version of the “field decay method” based on the cavity differential equation. The proposed approach 
enables the precise calibration of fo.s even under mismatch conditions. We tested the proposed approach on the 
SRF cavities of the Chinese Accelerator Driven System Front-End Demo Superconducting Linac and compared 
the results with those obtained from a network analyzer. The two sets of results were consistent, indicating the 
usefulness of the proposed approach. 
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I. INTRODUCTION 


Driven by the growing demand for safe nuclear fuel post- 
treatment processes, the China initiative Accelerator Driven 
System (CiADS) is being constructed as a clean solution for 
nuclear fission power sources [1-3]. To showcase the poten- 
tial of a high-power continuous wave (CW) proton beam for 
this project, the China ADS Front-End Demo Linac (CAFe) 
was built. This Linac is a 162.5 MHz superconducting (SC) 
radio-frequency (RF) machine operating in the CW mode 
and consists of both normal conducting (NC) and SC sec- 
tions (Fig. 1). The NC section includes an ion source, 
low-energy beam transport line, RF quadrupole accelerator, 
and medium-energy beam transport line. Conversely, the 
SC section comprises SC accelerating units, including 23 SC 
half-wave resonator cavities assembled into four cryomodules 
(CM1-CM4) [4-7]. The commissioning tests conducted on 
CAFe in the CW mode with a current of 10 mA and energy 
of 20 MeV successfully demonstrated its ability to accelerate 
and transmit high-intensity beams. 

For an SC cavity, the loaded quality factor (QL) reflects the 
consumption of the stored electromagnetic energy inside the 
cavity. In an ideal situation, in the absence of a beam pass- 
ing through the cavity, Q; includes the power dissipation on 
the cavity wall caused by the surface resistance (termed Qo 
[8, 9]) and the power leakage from the coupler ports [10]. 
Thus, Q; is a critical parameter that must be carefully se- 
lected to match the impedance of the RF generator with the 
particle beam load during operation [11]. Furthermore, dark 
current loading can negatively affect QL, making it an im- 
portant figure of merit for identifying such effects [12, 13]. 
In addition, Q; (or the cavity half-bandwidth (fo.5)) plays 
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a crucial role in the design of model-based controllers [14— 
16]. To satisfy the aforementioned application requirements, 
the measurement error for Qr, should not exceed 5%. The 
value of Q, can be calibrated using the cavity resonant fre- 
quency (fo) and fo.5, where fo.5 = ae Therefore, the pre- 
cise measurement of fo.5 is a prerequisite for calibrating QL. 
The decay curves of the cavity amplitude and phase (after the 
RF power is turned off) obtained from the cavity differential 
equation contain information on fo.5 and the cavity detuning 
parameter (A f), respectively. Many laboratories, including 
DESY, KEK, and CSNS, employ the “field decay method” to 
measure the aforementioned physical quantities [17—21]. 


We tracked the long-term changing regularity of fo.5 based 
on the data obtained when the RF power was turned off and 
accumulated while the CAFe facility was operated. Occa- 
sionally, we found that the cavity half-bandwidth calculated 
using the amplitude decay curve (i.e., fo.5,decay) and the 
cavity detuning parameter calculated using the phase decay 
curve (i.e., A faecay) appeared to be correlated. However, 
in principle, they should be independent. To better under- 
stand the above issue, we thoroughly examined the “field de- 
cay method.” Our findings revealed that this method is based 
on the zero-input response of the cavity differential equation, 
which indicates that the RF system must satisfy the “zero- 
input” condition. Thus, the cavity incident power must drop 
to zero after the RF power is turned off. However, if this con- 
dition is not met (i.e., owing to impedance mismatch), the re- 
maining incident power may influence the decay process and 
render the “field decay method” ineffective. We constructed 
an equivalent circuit that included RF power sources, trans- 
mission lines, input couplers, and SRF cavities and derived 
a solution for the cavity differential equation when a source 
impedance mismatch occurred. Finally, we modified the for- 
mula in the “field decay method” to explain the aforemen- 
tioned correlation. 


The “field decay method” is always employed to cali- 
brate Qg. Our findings revealed that this method is based 
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Fig. 1. Layout of the CAFe facility. Two types of half-wave resonator superconducting cavities (HWRO10 and HWR015) are implemented. 
The cavity CM3_3 is marked by a red triangle. Note that for a cavity CMm-n, the subscripts m and n represent the m*™™ cryomodule and 


nt” cavity, respectively. 


on the “zero-input” response of the cavity differential equa- 
tion. If the “zero-input” condition is not satisfied owing to an 
impedance mismatch, considerable errors may occur in the 
measurement of QL. To improve measurement accuracy, this 
study focuses on a modified calibration algorithm based on 
an equivalent circuit. 


Il. PHENOMENA AND POSSIBLE INTERPRETATION 
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Fig. 2. Schematic interpreting the notations of the cavity differential 
equation in polar coordinates. 


The conventional “field decay method” is briefly reviewed 
in this section. The naming rules for the cavity forward volt- 
age signal (Vç) and cavity voltage (Vo) in polar coordinates are 
illustrated in Fig. 2. In polar coordinates, Ve and V, can be 
expressed as Vp = pe? and V, = re’”, respectively, where 
p, r, 0, and represent the amplitudes of 2Vẹ and V. and 
the phases of V; and V., respectively. When the RF power is 
turned off, p immediately decreases to zero in the ideal case. 
Under this condition, according to the SC cavity differential 
equation without the beam in polar coordinates [22-24], fo.5 
and Af can be expressed as 


fo.5,decay =— (=) E : (1) 
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Fig. 3. Cavity amplitude (a and b) and phase (c) decay curves before 
and after the RF power is turned off for various Af on the cavity 
CM3-3. The parameter fo.5,decay is obtained by calculating the 
slope of the decay curve between 0.95 and 0.75 of the steady-state 
Ve (see the middle plot of c ). 


We refer to the various methods used in this study to cal- 
culate fo.5 and Af. For clarity, we first provide their detailed 
definitions in Table 1. 

For the cavity CM3_3 (marked by a red triangle in Fig. 
1) at CAFe, we measured the cavity-field decay curves for 
different values of A faccay (Fig. 3). fo was tuned using 
a frequency tuner. After the RF power was turned off, the 
slope of the cavity phase varied with the detuning parameter 
(Fig. 3(c)) because the phase decay curves are directly asso- 
ciated with the detuning parameter according to Eq. (1). Be- 
cause the cavity Qi is independent of the detuning parameter, 
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Fig. 4. (a) fo.5,decay and A faecay have dependency relations for different cavity voltages on CM3_3, particularly when Ve is greater than 0.4 
MV. (b) The dependency relations can also be observed in CM3~_4 (red triangles). (c) Comparison of the dependency relationships at CM3_3 
when A faecay is scanned by the signal source (red triangles) and tuner (gray dots). The two curves exhibit the same tendency. 


TABLE 1. Definitions of the cavity-related notations. 


Notation Definition 


fo.s 
fo.5, decay Cavity half bandwidth calculated using the field decay curve with Eq. (1) 


Actual cavity half bandwidth 


fo.5, cali Cavity half bandwidth calibrated using Eq. (15) 
0.5, scan Half bandwidth obtained by the network analyzer (scanning) 
Af Actual cavity detuning 
A foecay Cavity detuning calculated by the field decay curve with Eq. (1) 
A feai Cavity detuning calculated using Eq. (15) 
A fis Cavity detuning calculated from Ve and V+ in the steady-state 


the field decay curves of the cavity are expected to overlap 
under different detuning conditions; however, the they appear 
to be affected by the detuning parameter(Fig. 3(a) and (b). 
We conducted several studies to address this perplexing phe- 
nomenon. 

First, we calibrated fo.5 decay and A fdecay with four differ- 
ent values of V; for the cavity CM3_3 (Fig. 4(a)). The value 
of V, is less than the onset gradient of the field emission (ap- 
proximately 1.2 MV). All four curves show the dependencies 
between fo.5,decay and Afaecay. A similar dependence ap- 
pears in another cavity (CM3_,4) (Fig. 4(b)). 

Initially, we suspected that the frequency tuner might have 
disturbed the input coupler, causing variations in the coupling 
coefficients (8) and QL (or fo.5,decay). Consequently, we turn 
off the tuner and achieve cavity detuning by scanning the fre- 
quency of the signal generator at the same CM3_3. However, 
a similar dependence was observed in both cases (Fig. 4(c). 
The deviation in Fig. 4(c) is primarily because of the slight 
differences in the cavity field levels. 

Assuming that V;* and V,* represent the forward signal and 
reflected signal measured by the directional coupler, respec- 
tively, a further calibration of V;* and V,* is necessary to ob- 
tain the true forward and reflected signals (V; and V,) using 
Ve = XV," and V, = YV,*, respectively (assuming that we 


0.12 


Voltage [MV] 


Phase [deg] 


-100 1 f f 
0.5 1 1.5 


Time [ms] 


ML 


Fig. 5. Voltage and phase of Vc, XVř, and Y V,* before and after 
turning off the RF power. 


neglect the channel crosstalk). The complex coefficients (X 
and Y) can be obtained by solving a linear regression equa- 
tion [22]. Fig. 5 compares the V;, V; and V signals, and a 
residual attenuation signal of Vp (i.e., X V;") is observed after 
turning off the RF power. There are two possible reasons for 
this (Fig. 6). 


(1) The first reason is the crosstalk between the measure- 
ment channels, that is, the residual signal of V; is cou- 
pled with the signal V, [22]. In this case, the residual 
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Fig. 6. Two causes of V Æ 0 (RF off): crosstalk between the 
measurement channels (green line) and source impedance mismatch 
(red line). 


signal is a false measurement signal. 


(2) The second reason is the source impedance mismatch, 
where the V, signal is reflected from the generator side 
and mixed with the Vp signal. In this case, the residual 
signal is a true signal. 


Before 2021, directional couplers with poor directivity (ap- 
proximately 20 dB) were commonly used in our communica- 
tion system, which we refer to as CAFe. Previous studies 
have suggested that the limited directivity of these couplers 
was primarily responsible for the residual signals [24]. In 
2022, we replaced all the old directional couplers with new 
ones exhibiting high directivity (40 dB). This resulted in al- 
most negligible channel crosstalk; however, we decided not to 
install a high-power circulator in CAFe because of cost con- 
straints. Based on these factors, we conclude that the residual 
V; signal in Fig. 5 could be attributed to impedance mismatch 
rather than crosstalk. 

The algorithm in Eq. (1) must be modified because the 
“zero-input” condition is not satisfied. The specific calibra- 
tion algorithms are described in Sec. II. 


Il. THEORY AND ALGORITHM 


In this section, we establish cavity differential equations for 
the mismatched source impedance condition and use them to 
derive new formulas for calibrating fo.5 and Af. 


A. Radio-frequency and cavity circuit under the mismatched 
source impedance condition 


Fig. 7 (a) presents a simplified model of an RF cavity 
coupled to an RF generator using a rigid coaxial line and 
an RF input power coupler [19]. In this model, the coaxial 
line is represented by a transmission line with a characteristic 
impedance (Zo) and complex propagation constant (œ + iĝ). 
If the source impedance (Z,) (on the generator side) is not 


equal to Zp, a portion of the cavity-reflected signal is mea- 
sured by the direction coupler as the cavity forward signal 
after turning off the RF power. This process is described 
using an equivalent circuit (Fig. 7(b)). Assuming that the 
cavity input coupler has a transformation ratio of 1:N, the 
voltage signals V., Ve, and V,, can be transformed into V, 
Vi and vi. respectively, using the following transformation 
equation: V4 = + V. 

In Fig. 7 (b), the reflection coefficient (T`) for the forward 
generator side at z = 0 can be expressed as [25] 


V© Z- 
D= e 2) 
I V+ Z +Z 
where V, and V_ represent the voltages of the incident and 
reflected waves at z = 0, respectively. After the RF power 
is turned off, the reflection coefficient at z = —L (L is the 
length of the transmission line) is given by 


Vp  V-e-(0+iP)L 


_ —2(a+if)L 
L= UR > Wreariaye e0) 


T 


=Tye 


Therefore, after turning off the RF power, the transformed 
cavity voltage (VË) can be expressed as 


VE =v} + Vi, (4) 


Inserting Eq. (3) into Eq. (4) and eliminating vi. we obtain 


ly; 
vig ve. 5 
f 1 4 in Ë ( ) 
Accordingly, signal V¢ is associated with signal Vo by 
r 
V; = Vo. 6 
f 14+, (6) 


In the absence of the beam current, the cavity differential 
equation can be expressed as: 


dV. 
ae + (wo.5 — tAw)Ve = wo.5u, (7) 
and 
2be 
u = 2a Vp. (8) 


Here, the parameter 6e is the coupling factor, which is gener- 
ally considerably larger than 1 (i.e., 6, >> 1). Thus, u can 
be simplified as u = 2V;. By substituting Vf with Eq. (6) into 
Eq. (7), we obtain 


dV. ; 
We + (wo.5 — iAw)V: = wo.5s -a Vo, (9) 


where a = 2z /(1 + Tz) denotes a complex factor. Eq. (9) 
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Fig. 7. (a) Simplified model of a cavity coupled to an RF generator by a rigid coaxial line and an RF input power coupler. (b) Equivalent 
circuit diagram of (a) after the RF power is turned off. The cavity voltage (Ve) is transformed into V on the left side of the input coupler. 
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Fig. 8. Calibration error of the cavity half bandwidth (a) and detuning parameter (b) as a function of the reflection coefficient (TL), wherer, 
ranges from 0 — 0.18. The cavity bandwidth (fo.5,decay) and A faecay are calibrated based on the field decay curves. To ensure that the 
accuracies of fo.5,decay and A faecay / fo.5,caii lie within the +5% band, Fr should be located inside the red circle (i.e., || <0.024 ~ -32 
dB), as illustrated in (c). The results in (d) show that the calibration error of fo.5,decay and A faecay increases with lTxI, particularly when 
|['L| exceeds -32 dB. 


does not satisfy the “zero-input” condition. However, by re- following condition: 
arranging the terms, we obtain an equation that satisfies the 


dVe(t) 
dt 


+ [wo.s(1—a) —iAw(t)]Ve(t)=0, (10) 


B. New calibration algorithm for fo.; and Af 


It is more convenient to normalize the steady state (V.(t)) 
to one, that is, V.(0) = 1 if we assume that the RF power is 
turned off at t = 0 and the signal V(t) is in steady state at 
t < 0. Under the aforementioned restrictions, the solution to 
Eq. (10) at t > 0 is given by: 


V(t) = V.(0)eso [wo.5(a—1)+iAw(r)]dr 
(11) 


= alo [wo.5(a—1)+iAw(r)]dr f 


The complex factor (œa) can be decomposed into real and 
imaginary components (a, and aj), i.e., a = a; + 104. 

Separating Eq. (11) into real and imaginary parts yields the 
following: 


V(t) = e(or—1)wo.stt+i(aiwo.stt fo Aw(r)dr) (t > 0), (12) 
Consequently, the amplitude and phase of V; are given by 


T= |V.(t)| = e(ar—lwo.st 
(13) 


Eq. (13) provides some insights. If the source impedance 
(Zg) perfectly matches the transmission line impedance (Zo), 
then l'e = 0, T, = 0, and a = 0. In this scenario, accord- 
ing to the field decay algorithm, the cavity half bandwidth 
(wo.5) and cavity detuning parameter (Aw) can be easily ob- 
tained by fitting the slope of the cavity amplitude and phase, 
i.e., wo.5 = —r’/r and Aw = vy’. However, if the source 
impedance is not correctly matched to the transmission line 
(Tg Æ 0), the waves reflected from the cavity side will reflect 
again, resulting in nonzero forward power even after the RF 
power is turned off. This significantly affects the shape of the 
field curves, rendering the original algorithm inapplicable. 

Using the original field decay method, the cavity half- 
bandwidth (wo.5 decay) and detuning parameter (Awgecay) 
were calibrated using 


W0.5,decay = =” = (1 = Qr )wo.5 
(14) 


t 
Awdecay =P = Qiwo.5 + Aw 


To calibrate the actual cavity half-bandwidth and detuning pa- 
rameters, the original algorithm must be modified as follows: 


—r'/r _ fo.5,decay 
l-ar — l-ar 


fo.5,cali = 55° 

A feali — xy! = 04% fo.5,cali = A faecay = 0% f0.5,cali 
(15) 
The estimation error of the original field-decay algorithm 
can be easily evaluated for a specified L. According to 


Eq. (15), the accuracies of fo.5 decay and A faecay can be ex- 
pressed as: 


fo.5,decay—fo.5,cali __ 


2 207, 
Qr = r( 2) 


fo.5,cali 
(16) 
Afaecay-Afeai _ . _ a f 203 
fo.5,cali Bo led 14+P. 


The accuracy of the original field-decay algorithm depends 
on the parameter a, where the real (a) and imaginary (a;) 
parts of a determine the accuracies of fo.5 decay and A faecay; 
respectively. Fig. 8 (a) and (b) illustrate the calibration errors 
of fo.5,decay and A faccay, respectively, as functions of Tg. To 
ensure that the accuracies of fo 5,decay and A faecay / fo.5,cali 
lie within the +5% band, the coefficient r must be located 
inside the red circle (that is, |[,| <0.024 ~ -32 dB) in Fig. 8 
(c). Fig. 8 (d) illustrates the specific relationship between the 
accuracies of fo.5,decay> A fdecay / fo.5,cali and Tr 


IV. MODELING AND SIMULATION 


The state-space formalism in Eq. (7) is given by [19, 26, 
Z|: 


d Vor _ —W0.5 —Aw Ver 
dt Voi Aw —wo.5 Vai 
w 0 Ur 
4 [8 l (17) 
0 wos Uj 


Here, Ve r and Vo; represent the real and imaginary parts of 
the complex quantity (V.), while u, and u; represent the real 
and imaginary parts of the complex quantity (u), respectively. 
Eq. (17) can easily be transformed into its discrete-time form 
as follows [26, 27]: 


Vorn +1)| _ | 1-Tswos —TsAw(n +1)] | Ver(n) 
A +1) T,Aw(n +1) 1—T.wos Vea (n) 
T; 0 r 
WO0.5 u (n) s (18) 
0 T'sW0.5 ui(n 


where T, represents the sampling period. 
Eq. (9), the drive signal u can be expressed as 


According to 


i| |Ver(r 4+ 1) 
= (19) 
ui(n + 1) Qi Qr Vei(n +1) 


Other difference equations related to the dynamic behav- 
ior of the cavity are derived in Appendix A. The simulation 
parameters are presented in Table 2. After steady-state oper- 
ation is achieved for 0.4 ms, the RF power is turned off. The 
signal reflected from the cavity is also reflected as I’, is not 
zero. Fig. 9 compares the signals V., Ve, and V, based on the 
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Fig. 9. Comparison of the cavity voltage (a), cavity forward (b), and reflected signal (c) measurements based on the cavity model (red) and 
real SC cavity (gray). A stable operation is maintained for 0.4 ms before the RF power is turned off to create a field decay event. In the cavity 
model, the LFD dynamics is assumed to be determined by a first-order differential equation, resulting in curved trajectories in the cavity 
phase signal (red dashed lines) instead of linear ones. When the LFD effect is neglected, the cavity phase curve is linear (blue dotted lines), 
as shown in (a). In addition, the simulation results for the perfect impedance match case (Tẹ = 0, indicated by green) are also included for 
comparison. 


cavity model ( red) and real SC cavity (gray). Both results 


showed good consistency. In addition, the simulation results Te ee Ene ore oulsuons 


for the perfect impedance matching case (Tẹ = 0, indicated Item Value 

in green) are included for comparison. fz 184 Hz 
A finitial 149 Hz 
jal 0.31 
Za -42 deg 
Tm 530 us 


According to Appendix A, the LFD dynamics are assumed E wanes 
to be determined by a first-order differential equation. Con- al i 
sequently, the cavity phase signal presents curved trajectories 0.15 Hz/(MV/m)? 
(red dashed lines) rather than linear trajectories. However, Ts 80 ns 
without the LFD dynamics, the cavity phase curve is linear, 
as indicated by the blue dotted lines in Fig. 9(a). The slope of 
the linear phase curve represents the cavity detuning parame- 
ter (Awdecay). The cavity phase curves overlap in the first 80 
us after the RF power is turned off, regardless of whether the 
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V. EXPERIMENTAL VERIFICATION 


Based on the algorithm in (15), we recalibrated the mea- 


LFD is included (Fig. 9(a)). For clarity, the LFD dynamics 
during the field decay process were examined (Fig.10). For 
the subsequent 80 u s after the RF power is turned off, the 
LFD-induced phase error and maximum LFD value are less 
than 0.01 deg and approximately 0.8 Hz, respectively. Con- 
sequently, it is reasonable to fit the 80 u s cavity phase data 
to obtain Af. 


surement results shown in Fig. 4(a). The details of the process 
are described below. 


(1) Calibrating the actual forward and reflected signals us- 
ing Ve = X Vč and V, = XV,*, respectively: The com- 
plex factors (X and Y) were determined by solving a 
linear regression equation [22]. 
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Fig. 10. Simulation results showing the effect of the LFD on the 
measurement accuracy of the cavity phase and cavity detuning. For 
the initial 80 us period after turning off the RF power, the LFD- 
induced phase error is less than 0.01 deg, which can be ignored. 


(2) Determining the factor a: According to (6), we calcu- 
lated a as twice the ratio of V; to Ve after turning off 
the RF power. Here, we averaged a over a time inter- 
val of 50 us (e.g., from 0.45 ms to 0.5 ms in Fig. 11) to 
reduce uncertainty. 


(3) Calibrating fo.5 decay and A faecay using the traditional 
“field decay method”: To ensure a sufficient signal-to- 
noise ratio, we determined fo.5 decay by calculating the 
slope of the decay curve between 0.95 and 0.75 in the 
steady-state (V.)( Fig. 3(b)). To avoid the LED effect, 
we determined the derivative of the cavity phase within 
an 80 us interval after the RF power was turned off to 
obtain A fdecay- 


(4) Calibrating fo.s5cai and A feai: We used the for- 
mula (15) with the known values of a, fo.5,decay, and 
A fdecay to calibrate fo.5,cali and A fealis 


We used the aforementioned procedure to calibrate fo.5 cali 
and A fcaii (Fig. 12(a)). In contrast to the strong correlation 
observed between fo.5,decay and Afaccay measured by the 
cavity amplitude and phase decay curve, the new calibrated 
value of fo.5,ca Was not independent of A feaj at different 
Vo levels. 

We also compared A feaj and Afaccay with the steady- 
state cavity detuning parameters (A fss), which were calcu- 
lated from the steady-state values of Ve and V, (Fig. 12(b)). 


The following expression was used to calibrate A fis: 
Afss = tan(y — 9) - fo.5,cali- (20) 


Fig. 12(b)indicates that the discrepancy between A fecali and 
A fss is less than + 5 Hz, whereas there is an offset of ap- 
proximately 40 Hz between A faecay and A fies. 

To further validate the proposed algorithm, we used a net- 
work analyzer to measure CM3_3 and perform a compar- 
ison. Fig. 13(a) and (b)show the measurement setup, and 
the frequency-response measurement and corresponding fit- 
ted curve, respectively. The measurement results were fitted 
using the following formula 


Amax E k 
Aa = fos , (21) 


VF a Sottset } + fès 


where Amax is the maximum value of the amplitude- 
frequency-response curve. The parameters f and foffset 
represent the stimulus frequency and frequency offset be- 
tween the cavity resonant frequency and RF, respectively. In 
Fig. 13(b), the estimated value of fo.5 for the fitted curve is 
183.7 Hz. 

Next, we utilized a network analyzer to scan the cavities in 
CM, — CM; and estimated fo.5,scan by curve fitting. We then 
compared the derivation between fo.5 cai and fo.5 scan for the 
16 cavities (excluding CM>_¢ and CM3_, because of cavity 
faults). In addition, we plot the derivation between f.5 decay 
and fo.5 scan for comparison purposes. The results are shown 
in Fig. 14. The deviation between fo.5 scan and fo.5,cali Was 
maintained roughly within +2%, indicating that the proposed 
calibration algorithm accurately estimated the values of fo.5. 


VI. CONCLUSION 


After the RF power was turned off, the residual signal (V;) 
caused by source impedance mismatch could affect the field 
decay process, resulting in measurement errors for fo.5 and 
Af. The simulation results indicate that, to ensure that the 
measurement errors of fo, and Af were less than 5%, the 
reflection coefficient at the generator side must be less than 
-32 dB. To improve the measurement accuracy, we derived a 
calibration algorithm based on the cavity differential equation 
for cases in which there was a source impedance mismatch. 
Using this algorithm, we recalibrated fo.5 and Af and com- 
pared the results with those obtained from network analyzer 
measurements. The maximum error between the calculated 
and measured values was within 2%. 

q is a crucial factor in the calibration algorithm that ulti- 
mately depends on the impedance Z, on the generator side. 
Unfortunately, we do not have a comprehensive understand- 
ing of the factors that determine Z,. One possible factor is 
the power output of the RF generator. For instance, differ- 
ent V and detuning parameters require different RF powers, 
which lead to different 7, values. Nevertheless, it is impor- 
tant to stress that Z, is not determined solely by the generator 
power, and the same power may result in different values of 
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Fig. 11. (a) Cavity forward signal (V;) for different detuning values on cavity CM3_3(the corresponding cavity voltage signal is given in 
Fig. 3). (b) Method for estimating a from the Ve and V¢ signals after the RF power is turned off. 


Zg. In the next step, we focus on identifying the mechanisms 
and physical quantities that affect Z,. 
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Appendix A: Cavity dynamic behavior 


The LFD must be considered because the cavity gradi- 
ent varies after the RF power is turned off. As the cavity 
field attenuates slowly owing to the high loaded Q (narrow 
bandwidth) characteristic, higher-order mechanical modes 
are usually not excited. In this case, the LFD dynamics can 
be described by a first-order differential equation [19] as fol- 
lows: 


d 
Tm Afiep (t) + Aftro(t) = —Krro Epeak (AD 
where Ky,rp and 7 denote the LFD coefficient and me- 
chanical time constant, respectively. The quantity Epeak is 


the cavity peak gradient, which can be defined as V, using 
Epeak = Ve /Leg, where Leg is the effective cavity length. 


According to Ref. [27], the corresponding difference equa- 
tion of Eq. (A1) is given by 


Afurp(n +1) = bAfiep(n) + g [V2,(n) + V2,(n)] , 
(A2) 


where the coefficients are b = e~7*/T-FP and g = (1 — 
b), When the RF power is turned off, the cavity field 
ff 


exponentially decreases to zero. The corresponding LFD 
parameter value increases from its initial value of negative 
dozens ( hundreds) of Hz to zero. In addition, we assume that 
LFD is the only dominant cavity detuning that occurs after 
the RF power is turned off ( other sources may cause negli- 
gible detuning), and as V2 varies from V2(0) to V2(t), the 
LFD dynamics can be modified as follows: 


A furp(n +1) = bAfirp(n) +g [V2,(n) + V2;(n) - 1] , 
(A3) 


The constant | in Eq. (A3) represents the square of the 
steady-state (V.) amplitude, which is normalized to one for 
simplicity, that is, V?,(0) + V2,(0) = 1 

Thus, the overall cavity detuning is given by: 


Aw(n + 1) = 27 [Afinitia + Aftrp(n+1)].  (A4) 


Here, A finitia) is the initial steady-state detuning parameter 
before the RF power source is turned off. Using the difference 
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Fig. 12. (a) Relationship between the cavity half bandwidth (fo.5) and cavity detuning parameter (A fecali) in the calibrated and uncalibrated 
algorithms for various values of Vo. (b) Comparison of the results of A faecay and A f-aii measurements with the cavity detuning value (A fss) 
obtained in the steady state. 


equations (17), (19), (A3), and (A4), a dynamic simulation of 
the cavity undergoing a complete field decay process can be 
easily performed. 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


W.L. Zhan, Accelerator driven sustainable fission energy, in 
Proceedings of the 7th International Particle Accelerator Con- 
ference, IPAC2016, Busan, Korea, 2016. https://doi. 
org/10.18429/JACoW-IPAC2016-FRYAA03. 

Y.S. Qin, et al., Transfer line including vacuum differential 
system for a high-power windowless target, Phys. Rev. Accel. 
Beams 23 (2020) 113002, https://doi.org/10.1103/ 
PhysRevAccelBeams.23.113002. 

S.H. Liu, Z.J. Wang, H. Jia, et al., Physics design of the 
CIADS 25 MeV demo facility, Nucl. Instrum. Meth. A. 843 
(2017). https://doi.org/10.1016/4.nima.2016. 
10,055, 

Y. He, T. Tan, A.D. Wu, et al., Operation experience 
at CAFe, in Oral Presenta-tion of the 2021 International 
Conference on RF Superconductivity, SRF’21, in: vir- 
tual conference, JACoW, virtual conference, 2021. Avail- 
able athttps://indico.frib.msu.edu/event/38/ 
attachments/160/1298/MOOFAV03_yuan_he.pdf 
Q. Chen, Z. Gao, Z. L. Zhu, et al., Multi-frequency point sup- 
ported LLRF front-end for CiADS wide-bandwidth applica- 
tion, Nucl. Sci. Tech. 31 (29) (2020). https://doi.org/ 
10. 1007/5s41365-020-0733-9 

W.M. Yue, S. X. Zhang, C. L. Li, et al., Design, fabrication and 
test of a taper-type half-wave superconducting cavity with the 
optimal beta of 0.15 at IMP, Nuclear Engineering and Technol- 


[10 


(ll 


[12 


] 


ol 


REFERENCES 


ogy, 52 (8) (2020): 1777-1783., https://doi.org/10. 
1016/j.net.2020.01.014. 

W.M. Yue, S.X. Zhang, C.L. Li, et al., Development of a low 
beta half-wave superconducting cavity and its improvement 
from mechanical point of view, Nucl. Instrum. Meth. A. 953 
(2020) 163259, https://doi.org/10.1016/j.nima. 
2019.163259; 

P. Sha, W.M. Pan, S. Jin, et al. Ultrahigh accelerating gradient 
and quality factor of CEPC 650 MHz superconducting radio- 
frequency cavity, Nucl. Sci. Tech. 33(10):125(2022).ht tps : 
ffdol.erg/ 10.1007 /s41365-022-01109-8 

Z.Y. Ma, S.J. Zhao, X.M. Liu, et al. High RF power tests of 
the first 1.3 GHz fundamental power coupler prototypes for the 
SHINE project. Nucl. Sci. Tech. 33(1):10(2022). https: // 
doi.org/10.1007/s413¢65-022-00964=5 

J.P. Holzbauer, C. Contreras, Y. Pischalnikov, et al., Improved 
RF measurements of SRF cavity quality factors, Nucl. Instrum. 
Meth. A. 913 (2019) 7-14 https://doi.org/10.1016/ 
Jsnima,2018.09.155 

F. Marhauser, Method for in situ and in operando cavity loaded 
Q extraction in superconducting rf accelerators, Phys. Rev. Ac- 
cel. Beams, 24, 032001 (2021), https://doi.org/10. 
1103/PhysRevAccelBeams.24.032001. 

A. Bellandi, L. Butkowski, B.Dursun, et al., Online Detuning 
Computation and Quench Detection for Superconducting Res- 


Cavity 


i 
LN 


J 


Network analyzer 


(a) 


Mag. [dB] 


Phase [deg] 


11 


50 


oO 


a 
ro) 


100 L 1 L i 1 L 
-500 -400 -300 -200 -100 0 100 
Frequency [Hz] 


200 300 400 500 


(b) 


Fig. 13. (a) Measurement block diagram based on the network analyzer. (b) Frequency response function of cavity CM3~3, as measured by 
the network analyzer. The foffset parameter is excluded for better comparison. 


4 

— ob HE with calibration: (fp 5 ai40.5,scan!fo.5,scan 

x 

5 0 

Wi oL O E J 
-4 
40 T T T T T T 
30 > |E wo calibration: (o.s.decayf0.5,scan!fo.5,scan | 
20 | HE With calibration: (fp 5 cai fo.5,scan)fo.5,scan | 

L 

o 

5 10 F L 4 

i 

E T 

2 

= | 

< -10 | [ f [ | | | 

o 

i -20 + J 
-30 + 4 
-40 bie er, ee ge ge ej 


Me Bs ee I eR ae a 


1 1 1 
[oe] + [to] o 
Oo oo 9 9 
=>asaease 
O O O O 


CM, 
CM, 
CM,.3 
CM, 4 
M 
M 
M 
M 
M 
M 
M 
M 


Cavity Number 


Fig. 14. Comparison of the deviation between fo.5 decay and fo.5,cali 
for the 16 cavities and the network analyzer measurement results, 


fo.5,cali- 


onators, IEEE Transactions on Nuclear Science, vol. 68, no. 4, 
pp. 385-393, April 2021, https: //doi.org/10.1109/ 
TNS .2021.3067598. 

[13] J.Branlard, V.Ayvazyan, O.Hensler, et al., Superconduct- 
ing Cavity Quench Detection and Prevention for the Euro- 


pean XFEL, in Proceedings of ICALEPCS2013, San Fran- 
cisco, CA, USA. Available at https://accelconf.web. 
cern.ch/ICALEPCS2013/papers/thhppc072.pdf 

[14] F. Qiu, S. Michizono, T. Miura, et al., Application of distur- 
bance observer-based control in low-level radio-frequency sys- 
tem in a compact energy recovery linac at KEK, Phys. Rev. 
ST Accel. Beams. 092801 (2015), http://doi.org/10. 
1103/PhysRevSTAB.18.092801 

[15] F. Qiu, S. Michizonol, T. Matsumotol, et al., Combined 
disturbance-observer-based control and iterative learning con- 
trol design for pulsed superconducting radio frequency cav- 
ities, Nucl. Sci. Tech. (2021) 32:56, https://doi.org/ 
10.1007/s41365-021-00894-y 

[16] Y. Cong, S.F. Xu, W.X. Zhou, et al., Low-level radio-frequency 
system upgrade for the IMP Heavy Ion Reach Facility in 
Lanzhou (HIRFL), Nucl. Instrum. Meth. A. 925(2019):76- 
86. https://doi.org/10.1016/j.nima.2019.01. 
086. 

[17] M. Omet, Digital Low Level RF Control Techniques and Pro- 
cedures Towards the International Linear Collider (Ph.D. the- 
sis), School of High Energy Accelerator Science, 2014. 

[18] F. Qiu, T. Miura, D. Arakawa, et al., RF commissioning of 
the compact energy recovery linac superconducting cavities 
in pulse mode, Nucl. Instrum. Meth. A. 985 (2021) 164660, 
http://doi.org/10.1016/j.nima.2020.164660. 

[19] T. Schilcher, Ph.D. thesis, Universitt Hamburg, 1998. 

[20] J.Y. Ma, G.R. Huang, Z. Gao et al., The resonant frequency 
measurement method for superconducting cavity with Lorentz 
force detuning. Nucl. Instrum. Methods A 993(5), 165085 
(2021). https://doi.org/10.1016/j.nima.2021. 
165085 

[21] Z.C. Mu, Overview of the CSNS linac LLRF and operational 
experiences during beam commissioning, in Proceedings of 
the 57th ICFA Advanced Beam Dynamics Workshop on High- 


[22] 


[23] 


[24] 


Intensity, High Brightness and High Power Hadron Beams, 
Malmo, Sweden, 2016, pp. 3-8. 

A. Brandt, Development of a Finite State Machine for the Au- 
tomated Operation of the LLRF Control at FLASH(Ph.D. the- 
sis), University Hamburg, 2007. 

M. Grecki, V. Ayvazyan, J. Branlard, et al., On-line RF am- 
plitude and phase calibration for vector sum control, in Oral 
Presentation of the 9th Low-Level RF Workshop, LLRF2019, 
Chicago, USA September 29 — October 3, 2019. 

F. Qiu, J.Y. Ma, G.D Jiang, et al., Approach to calibrate ac- 
tual cavity forward and reflected signals for continuous wave- 
operated cavities, Nucl. Instrum. Methods A 1034, 166769, 


[25] 
[26] 


[27] 


12 


(2022). https://doi.org/10.1016/j.nima.2022. 
166769. 

David, M. Pozar., Microwave Engineering, 4th Edition, Wiley. 
T. Czarski, K.T. Pozniak, R.S. Romaniuk, and S. Simrock, 
TESLA cavity modeling and digital implementation in FPGA 
technology for control system development, Nucl. Instrum. 
Methods Phys. Res., Sect. A 548, 283 (2005). https: // 
doi,o6rg/10.1016/7.rima.2005.10.122 

F. Qiu, S. Michizono, T. Miura et al., Real-time cavity simu- 
lator based low-level radio frequency test bench and applica- 
tions for accelerators, Phys. Rev. Accel. Beams 21, 032003 
(2018). https://doi.org/10.1103/PhysRevSTAB. 
18.0928014 


